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! ABSTRACT: For the embedded Gaussian orthogonal ensemble (EGOE) of random ma- 
trices, the strength sums generated by a transition operator acting on an eigenstate vary 

T-H ■ with the excitation energy as the ratio of two Gaussians. This general result is compared 

, to exact shell model calculations, with realistic interactions, of spherical orbit occupancies 

\^ [ and Gamow- Teller strength sums in some (ds) and (fp) shell examples. In order to con- 

■ firm that EGOE operates in the chaotic domain of the shell model spectrum, calculations 
. are carried out using two different interpolating hamiltonians generating order-chaos tran- 
I sitions. Good agreement is obtained in the chaotic domain of the spectrum, and strong 
' deviations are observed as nuclear motion approaches a regular regime (transition strength 
• sums appear to follow the Dyson's A3 statistic). More importantly, they shed new light 

I , on the newly emerging understanding that in the chaotic domain of isolated finite interact- 

CJ ' ing many particle systems smoothed densities (they include strength functions) define the 

2 ■ statistical description of these systems and these densities follow from embedded random 

. . , matrix ensembles; some EGOE calculations to this end are presented. 

ixl ■ 

■ 1. Introduction 

. P.: 

In the last fifteen years there has been an explosive growth in the use of random matrix 
theories for quantum systems particularly in the context of quntum chaos |]l|, ^ . Recently, 
working with the aim of developing a statistical theory for finite interacting many parti- 
cle systems, such as atoms, molecules, nuclei, atomic clusters, metallic quantum dots etc., 
by incorporating the ideas of random matrices and chaos, several research groups recog- 
nized the importance of investigating the embedded random matrix ensembles in detail, 
i.e. the embedded Gaussian orthogonal ensemble of random matrices of /c-body interac- 
tions (EGOE(A;)) and their various deformations [2-7]. The EGOE is introduced in the 
context of nuclear shell model studies |2[; The EGOE(/c) is defined in m-particle spaces 
(i.e. in the ^ ) dimensional space generated by distributing m fermions over Af single 

particle states) with a GOE representation in A:-particle space for A:-body operators (usually 
k « m). This ensemble and its deformed versions are studied using Monte-Carlo methods 

^Invited talk in the National Seminars on Nuclear Physics held at Institute of Physics, 
Bhubaneswar during July 26-29, 1999. 
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with EGOE constructed in occupation number representation and alternatively assuming 
that realistic nuclear (similarly atomic) shell model hamiltonians are typical members of 
EG0E(2) and its deformations, generic features of these ensembles are inferred from exact 
shell model calculations. For finite interacting particle systems in the chaotic domain there 
are now several studies of the statistical properties (both smoothed forms and fluctuations) 
of energy levels, wavefunction amplitudes or equivalently transition strengths generated by 
action of a transition operator on an eigenstate and strength functions [1-11]- However, only 
in the last 2-3 years similar studies (i.e. in the context of quantum chaos) on expectation 
values of operators, which measure transition strength sums, as function of excitation en- 
ergy have began [3,8,12-14]. Given an operator K = O^O, the expectation values (K)^ are 
the diagonal elements of K in the hamiltonian (H) diagonal basis (a more precise definition 
is given ahead in Eq. (3)); they give strength sums generated by the transition operator O 
acting on the eigenstate with energy E. Two major examples are single particle transfer 
strength sums which are expectation values of number operators that give occupancies of 
single particle states (they determine thermodynamic behaviour) and Gamow- Teller (GT) 
strength sums as function of excitation energy which are relevant in astrophysics (presu- 
pernova evolution and stellar collapse). It is expected that the smoothed (K)^ vs E will 
give information about order-chaos transitions just as energies, wavefunction amplitudes 
and transition strengths. 

Two important results given by EGOE are that in strongly interacting shell model spaces 
(essentially in Ohuj spaces) (i) the state densities take Gaussian form and (ii) the bivariate 
strength densities take bivariate Gaussian form jl^. These results have their basis in the 
EGOE representation of the hamiltonian H (which is in general one plus two-body in nuclear 
case) and transition operators O. The eigenvalue density I{E) or its normalized version 
p{E) is defined by g 

I{E) = mH - E))) = d {5{H -E))=d p{E) ; 

EGOE^^ 1 lfE-e\'^ (1) 

p{E) — > p{E) = pg{E) = -j=^exp - - 



In (1) ((•••)) denotes trace (similary (• • •) denotes average), the e, o and d are centroid, 
width (fj^ is variance) and dimensionality respectively. Note that e = o"^ = ((i/ — e)^), 
stands for Gaussian and the bar over p{E) indicates ensemble average (smoothed) with 
respect to EGOE. The strength i?(-E, -E') generated by a transition operator O in the K- 
diagonal basis is i?(£', £") = | (£" \ 0\ E) p. Correspondingly the bivariate strength density 
hiv]oiE, E') or phiv-oiE, E') which is positive definite and normalized to unity is 



h^.,o{E,E') = (^(0U{H - E')05{H - E) 
= r{E') \{E'\0\E) |2 I{E) 



EGOE 



Pbiv;o{E,E') — > pbiv-o{E, E') — phi^_g.a{E, E') 
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The bivariate Gaussian phi^-g-oiE, E') in (2) is defined by the centroids and variances of its 
marginal densities and a bivariate correlation coefficient. Though the EGOE forms in (1,2) 
are derived by evaluating the averages over fixed-m spaces, in large number of numerical 
shell model examples it is verified that [2,10,15,16] they apply equally well in fixed- m, mT 
and mJT spaces. In practice the so-called Edgeworth corrections are added to the Gaussian 
forms in (1,2). Before proceeding further two remarks are in order: 

1. Level and strength fluctuations for EGOE follow the GOE fluctuations, i.e. nearest 
spacing distribution is Wigner distribution, Dyson-Mehta A3(L) statistic follow the 
GOE ln{L) behaviour for large L and strength fluctuations are of Porter-Thomas 
(P-T) type. The A3 form is tested ahead in Fig. 3 and recent tests of P-T form for 
EGOE are given in Fig. 1. 

2. EGOE smoothed forms (1,2) (and (3) ahead) gave birth to the so called statistical 
nuclear spectroscopy (see [9,10,15-17] and references therein) and there are recent 
studies of this in atoms |l^, molecules and solids |l^ and mesoscopic systems [3-5]. 



The pupose of this paper is to first point out that EGOE, via (1,2) gives rise to a statistical 
theory for the smoothed forms for transition strength sums and the theory operates in the 
chaotic domain of the spectrum. Shell model tests in [ds)^ space for occupancies and 
GT strength sums in {fp)^ space are carried out. The EGOE theory and the shell model 
studies are described in Sect. 2. In order to confirm that the agreement between shell 
model and EGOE theory is a consequence of chaoticity of the shell model spectrum, GT 
strength sums and occupancies are calculated using two different interpolating hamiltonian 
that generate order-chaos transitions. Results of these calculations form Sect. 3. Results 
in Sects. 2, 3 give nuclear physics examples for the newly emerging understanding that in 
the chaotic domain of isolated finite interacting many particle systems smoothed densities, 
defined by EGOE, give rise to the statistical description of these systems. Further remarks 
on this important generic result for quantum chaos in isolated finite interacting particle 
systems are given (together with the results from a EGOE calculation for occupancies) in 
the concluding Section 4. 

2. EGOE theory for transition strength sums and shell model tests 



One important by product of (2) is that the transition strength sum density {^{0^0 5{H - E) 
which is a marginal density of the bivariate strength density, takes a Gaussian form as the 
marginal of a bivariate Gaussian is a Gaussian. Therefore, using (1), it is immediately 
seen that transition strength sums vary with excitation energy as ratio of Gaussians. Given 
K = 0^0^ the transition strength sum density is the expectation value density denoted by 
Pk{E) and then ^ |l|], 

{Kf = [d{m)p{E)]-^ ^ (Sa I I Ea) = Ik{E)/I{E) = pk{E)/p{E) 

VaeE J 

^^^J^)/-^)=PK:g{E)/pg{E) , 
Pk{E) = {K6{H - E)) = d-^ Ik{E) = d'^ {{K6{H - E))) ; K = O^O . 



3 




X = Ft(Ei,Ef) / FC<Ei,Ef) 




Fig. 1. (a) Distribution of renormalized shell model transition strengths compared with 
the Porter-Thomas (P-T) form. The locally averaged strengths R{Ei, Ef) are calculated 
using the EGOE bivariate Gaussian form (2). Shell model (SM) calculations are in 223 
dimensional (cis)"^=5'"^=5/2.^=i/2 space with hamiltonian defined by Wildenthal interaction 
1^] and the transition operator is the two-body part of the hamiltonian obtained after 
substracting the configuration isospin centroid part. The SM+EGOE result is in good 
agreement with P-T. (b) Number of principal components (NPC) and information entropy 
^^m/o-j versus energy {E) for a strength distribution in 307 dimensional ((is)™=6,-'=2,T=o 
space. The hamiltonian is Kuo interaction with ^^O single particle energies ||2^ and the 
transition operator is same as in (a) but for the Kuo interaction. The exact shell model 
results are compared with the GOE and EGOE predictions; the EGOE formulas (they use 
P-T) are given in Q. (c) Same as (b) but for wavefunctions. 
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In deriving (3) it is assumed that the smoothed form of pk{E) / p{E) reduces to ratio of 
smoothed forms of Pk{E) and p{E). This result ignores the fluctuations in both pk{E) and 
p{E) and the r.m.s error due to neglect of the fluctuations is given in terms of the number 
of principal components (NPC) or the inverse participation ratio for the transition operator 
O [||]. Note that (3) takes into account {K,H) and {K.H'^) correlat ions which define the 
centroid eK and width gk of Pk{E); ex =< KH > / < K > and ctIk =< KH"^ > / < 



K > —e\. First discussions of the EGOE result in (3) are in [16, 17, 19]. Statistical models 
that are inappropriate, are applied recently in the study of GT strength sums as function 
of excitation enegy in (ds)-shell [ |T^ ] although there are several studies of GT strengths and 



strength sums in statistical nuclear spectroscopy pq, 21, 22 



In order to study the domain of validity of (3), for the occupancies (ria)^ shell model calcu- 
lations in 307 dimensional (2slci)'"^^''^^^'^^'^ space are carried out using the Rochester - Oak 
Ridge shell model code and for the GT strength sums in 814 dimensional (^lf2p)"^^^''^^^''^^^ 
space using the NATHAN code ||2^ of the Strasbourg-Madrid group. In the (ds)-shell stud- 



ies the hamiltonian employed is H = h{l) + V{2) defined by Kuo's |24] two-body matrix 



elements {V{2)) and ^''O single particle energies {h{l) <^=^> e^g^^ = —4.15 MeV, e^g^^ = 0-93 
MeV, £^^^2 = —3.28 MeV). Similarly in the (fp) shell studies the so called KB3 interaction 
H is employed with h{l) e/^^^ = 0.0 MeV, e/^^^ = 6-5 MeV, ep^^^ = 2 MeV, Cp^^^ = 4 
MeV). The expectation value density Pn^-.g for the number operators Ua in (ds) shell and 
PK(GT):g for the K{GT) operator that generates GT strength sums in {fp) shell are con- 
structed in terms of their centriods and widths and similarly the state density Gaussian. 
Then using (3) the smoothed form of GT strength sum as function of excitation energy 
is constructed and compared with exact shell model results. From Fig. 2 it is seen that 
the EGOE result (3) describes very well the shell model results except at the edges of the 
spectra. The reason for the deviation at the edges is well known - here the states are not 
chaotic (sufficiently complex). In the {ds) shell example the X-density centroid ek, width 
uki skewness ^i-k and excess '^2:K (the 71 and 72 are shape parameters that measure de- 
viations from the Gaussian form) are ex = — 35.08MeV, ax = 9.63MeV, '^i-k = 0.08 and 
72:j^ = -0.09 for ^5/2 density and (-29.50 MeV, 10.24 MeV, -0.08, -0.17) respectively 
for the (i3/2 density. Similarly for the state density the parameters are e = — 32.78MeV, 
a = 10.24MeV, 71 = 0.05 and 72 = -0.18. In the {fp) shell example, ex = 11.12MeV, 
ax = 8.65MeV, ^h-.k = 0.09, 72:^ = -0.18, e = 9.51MeV, a = 8.62MeV, 71 = 0.10 and 
72 = —0.19. Firstly the |7i| and I72I values (being much less than 0.3) clearly show that 
all the densities are close to Gaussian. Moreover a = ax/f^ ~ 1- The scaled centroid shifts 
= i^K — e)/c are —0.225 and 0.32 for the d^/2 and ^3/2 densities while it is 0.187 for the 
{fp) shell GT example. With {K)^ given as ratio of Gausians and that a ~ 1 imply that 
in the middle of the spectrum (K)^ ~ (K) |l + ^kE^', E = {E — e)/a, i.e (K)^ is linear 



in energy. This linear form and its polynomial extensions are used in the past |2^, 25 
The value of E for the ground state are —2.6 and —3.1 for {ds) and {fp) shell examples 
respectively. The range of validity (from the center of the spectrum) of the linear form for 
{K)^ is much smaller for ((is)-shell than in (/p)-shell as for the former is larger than 



for the later; see [22|. Fig. 2 showing good agreement between shell model results for (n^)^ 



and {K{GT))^ with realistic interactions and the EGOE prediction (3), gives a justifica- 
tion to our assertion that EGOE starts operating from the region of the onset of chaos for 
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the strength sums as well. Next we study this question via (ua)^ and {K{GT))^ when 
the hamiltonian changes thorough a parameter from a symmetry preserving, i.e. regular 
hamiltonian to a chaotic one. 




-60 -40 -20 
(a) Energy (MeV) 



Fig. 2. (a) Occupation numbers for ^5/2 and ^3/2 orbits vs excitation energy {E) in the 
same example as in Fig. lb. Shown are also the shell model state and occupation number 
densities (histograms) compared to the EGOE Gaussian forms (continuous curves) given 
by (1,3). (b) GT strength sum versus excitation energy (E) for the 814 dimensional six 
particle (/p)-shell space with J = 0, T = 0. The exact shell-model results for the realistic 
KB3 interaction |23| are compared with the EGOE predictions given by (3). 



3. Results v^ith MF and SU(4) interpolating hamiltonians 
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For further confirming the conclusions from Fig. 2, shell model calculations are performed 
in the 325 dimensional (ds)'^^^''^^'^^^ space for occupancies and GT strength sums with 
two different interpolating hamiltonians . First set of calculations use the spherical shell 
model mean-field (MF) hamiltonian h{l) as the unperturbed hamiltonian Hq and in this 
case the occupation number operators commute with Hq, 



Hx{MF) = hil) + Ay (2) = Ho + \{Hsm - Hq) 



(4) 



Note that Hsm = h{l) + V{2); the is defined by ""^"^O single particle energies and V{2) 
by Kuo's two-body matrix elements as in the Sect. 2 ((is)-shell example. In the figures in 
Figs. 3a-d the calculations with (4) are denoted by MF. It is easily seen that spherical 
configurations (generated by distributing the nucleons in the three (ds) shell orbits) are 
eigenstates for Hq. Therefore for A = in (4), the spectrum will have degeneracies. In the 
second set of calculations the S'C/(4) — ST scalar part Hsi;{i)-ST of Hsm is used as the 
unperturbed hamiltonian, 



Hx{SU{^)) = H, 



SU{4:)-ST:scalar 



SM 



H 



SU (A) -ST -.scalar) 



Ho + X{HsM - Ho) (5) 



In the figures in Figs. 3a-d calculations with (5) are denoted by SU(4). Note that the 
GT operator commutes with the SU(4) hamiltonian -ffs(7(4)-5T- For H = i?sJ7(4)-ST) the 
eigenvalues and eigenvectors are given easily by 5'C/(4) — ST algebra. The eigenstates are 
labelled, for a given number m of valence nucleons, by L, S, J, T and the C/(4) irreducible 
representations (irreps) {/} = {/1/2/3/4} or the S'C/(4) irreps {F} = {F1F2F3} where 

/l + /2 + /3 + /4 = m, /l > /2 > /3 > /4 > 0, Fi = (/i + f2-fs- /4)/2, F^ = 

(/i - /2 + /s - /4)/2, F3 = (/i - /2 - /a + fi)/2 with additional restrictions on /j's coming 
from the spatial part. In addition, as GT operator is a generator of the 5^7(4) group, the 
SU (4) algebra gives directly {K{GT))^ in terms of SU (A) — ST quantum numbers. Methods 
for constructing Hgjj(^^^_gj. part of a given H are given in |2G], 



Hsu{A)-ST = i(n-l)(n-2)i^(0,{0}00)-n(n-2)i?(l,{l} 



1 i> 

2 2^ 



+- 



n(2n - 9) + G2 + 2(52 + r2)J F(2,{2}11) 
n - G2 + 2(5^ + T^)] F(2,{2}00) 



n(n + 3) - G2 + 2(5^ - T^) F(2,{11}10) 



(6) 



+ - n(n + 3) - G2 - 2(52 - T^) F(2, {lljOl) 
o L J 

In (6) n is number operator, G2 is U{A) quadratic Casimir operator with eigenvalues 
(^2)^^^ = + 3) + /2(/2 + 1) + /3(/3 - 1) + /4(/4 - 3) and 52 and are operators with 
eigenvalues S'(S'+1) and T(T+1). Construction of -f^5[/(4)-5T requires the values for the cen- 
triod energies (determined by the given H) E{m,{f}ST) = (H)""'^^^^^. For the ((is)-shen 
Hsm hamiltonian (defined after (4)), they are E{0,{0}00) = 0, E{1,{1}^^) = -2.302MeV, 
£;(2,{2}11) = -4.176MeV, £;(2,{2}00) = -2.975MeV, £;(2,{11}10) = -8.360MeV and 
-E(2,{11}01) = — 7.048MeV. Using these, H\ is constructed and then occupancies and 
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{K{GT))^ are calculated for various A values. The SU{4) — ST algebra gives for A = 0, 
(i^(GT))^ = |((C2(5C/(4)))^^>-5(5 + l)); (C2(SC/(4)))^^> = Fi(Fi + 4) + F2(F2 + 2) + F|. 
In addition, the allowed {f}S values for the ((is)-shell example are: {2222}0; {3221}!; 
{3311}2,0; {332}1; {4211}!; {422}2,0; {431}1,2,3; {44}0,2,4; {5111}0; {521}!, 2; {53}1,3; 
{611}!; {62}2. With these results, the GT curve for A = case is constructed as a check 
of the shell model calculations. Each of the {f}S states will have several L = S states and 
therefore the states here have degeneracies. 

It is clearly seen from the results in Figs. 3a, b that the EGOE smooth form is not a 
good approximation to the exact results in the case of regular motion. For A ~ there 
are several (approximately) good quantum numbers with nearby levels carrying different 
sets of quantum numbers and therefore expectation values show large fluctuations as a 
function of excitation energy. The order-chaos transition as A increases is clearly illustrated 
by the spectral regidity A3 in Figs. 3c, d (also by the distribution of nearest neighbour 
level spacings as shown in [14]). It is seen that occupancies and GT strength sums behave 
rather like the A3 statistic, approching more slowly the EGOE and GOE limits respectively. 
This similarity is probably due to the fact that both A3 and strength sums are related to 
long-range correlations between the energy levels or wavefunctions. Thus, in the quantum 
chaotic domain transition strength sums, independent of the hamiltonian, follow EGOE 
forms and we have statistical spectroscopy in the chaotic domain. 

There is another important observation that follows from Figs. 3a-d, i.e. as the interact- 
ing particle system becomes chaotic, expectation values take smoothed forms (within GOE 
fluctuations |Q) and hence described by the smoothed densities. For A >> Ac (Ac corre- 
sponds to order-chaos border and we determine Ac by using xfx) which is the mean square 
deviation of the exact shell model strength sum from the EGOE smoothed form [0]; see 
Fig. 4) they take the EGOE form given by (3). This generic result is of central interest in 
quantum chaos studies of finite interacting particle systems as discussed in the next section. 



4. Further results and concluding remarks 



Besides the nuclear shell model results for GT strength sums and occupation numbers 
presented in Figs. 2-4 (GT results for the SU(4) case are reported in [14] recently), the 
behaviour of occupancies as a many particle hamiltonian makes order-chaos transitions is 
studied recently by several groups: (i) using a 20 member EG0E(l-|-2) in 330 dimensional 
= ll,m = 4 space with the MF hamiltonian (4); /i(l) is defined by the single particle 
enegies ei = i + (l/i) ; i = 1,2, ...,11 and V{2) is EG0E(2) S; (ii) using the four 



interacting electrons Ce atom [|ll|; (iii) using a symmetrized coupled two-rotor model |p~3| ; 
(iv) using EG0E(l+2) as in (i) but in the 924 dimensional M = 12, m = 6 space with 25 
members [|2^. Most significant conclusion of all these studies is that transition strength 
sums show quite different behaviour in regular and chaotic domains of the spectrum. In 
order to make this argument clear results of (iv) are shown in Fig. 5; here occupancies for 
the single particle states are calculated for various values of the interpolating parameter A. 
It is clearly seen, for example from Fig. 5, that below the region of onset of chaos 
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Fig. 3. (a) Occupation numbers for ^5/2 and ^3/2 orbits vs excitation energy (E) for 
the MF and SU(4) interpolating hamiltonians given by (4) and (5). Shell model results are 
compared with the EGOE predictions (continuous curves) given by (3). (b) same as (a) but 
for GT strength sums, (c) Averaged spectral rigidity A^lL) for the eigenvalues of the MF 
hamiltonian (4). Error bars give the standard deviation of the A3 average over overlaping 
intervals of length L. The dashed curves are for Poisson and the continuous curves are for 
GOE. (d) same as (c) but for the SU(4) hamiltonian (5). For the A = 1 case the parameters 
(e, £7,71,72) for the state, c?3/2 and GT densities are (—52.59 MeV, 13.15 MeV, 0.10, 
0.03), (-55.30 MeV, 12.31 MeV, 0.01, -0.06), (-48.70 MeV, 13.42 MeV, 0.07, 0.05) and 
(-48.37 MeV, 12.64 MeV, 0.06, -0.1) respectively. 
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GT.GT 
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Fig. 4. R{\) = xfx)/'^'{i) Results for MF and SU(4) calculations are shown for 



^5/2 occupancies and GT strengtii sums. From the observations in the studies in a 
plausiable definition that the Hx system is chaotic is given by the condition i?(A) < 1.5. 
Then, Ac is defined by R{Xc) = 1-5. The Ac values for the four cases are shown in the figure. 
The difi'erence in the two Ac values in the case of occupancies is easily understood in terms 
of the norms |15, ^] of Hq and Hi = Hsm — Hq; however this is not the case with GT 
strength sums. Thus it appears that a complete theory for Ac, in case of strength sums, 
may come from the study of A3, strength fluctuations, NPC and information entropy for 
the hamiltonian and the transition operator involved. 
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EGOE(1+2) : (N=12,m=6), Hx^= h + :^V} 




12 -2-10 
(E-8)/a 



2 3 



Fig. 5. Occupation numbers for a 25 member EE(l+2) ensemble, defined by the hamil- 
tonian h{l) + X{V{2)}, in the 924 dimensional = 12, m = 6 space (see text); in the 
figure M is denoted as A^. Details of matrix construction etc. are given in ||2^. Results 
are shown for the lowest 5 single particle states and for six values of A. In the calculations 
occupation numbers are averaged over a bin size of 0.1 in (E — e)/a; e is centroid and a is 
width. The spectra of all the ensemble members are first zero centered and scaled to unit 
width and then the ensemble average is carried out. The estimate of Q gives Ac ~ 0.05 for 
order-chaos border in the present EE(l-|-2) example. It is clearly seen that once chaos sets 
in, the occupation numbers take stable smoothed forms. See |27] for further details. 
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occupation numbers show strong fluctuations (in the regular ground state domain pertur- 
bation theory applies). In this region there is no equilibrium distribution for IxiE), I{E) 
and other densities. However in the chaotic domain (see ^] and also Fig. 4 for methods 
of determining Ac) the densities can be replaced by their smoothed forms as in (l)-(3). 
Therefore there is a statistical mechanics, defined by various smoothed densities, operat- 
ing in the quantum chaotic domain of finite interacting particle systems (this is also the 



essence of statistical nuclear spectroscopy [9,10,15-17]; see |28] for arguments for statistical 
spectroscopy in atoms). In fact in favourable situations, it is possible to introduce ther- 
modynamic concepts (effective temperatures and chemical potencials etc.) in the chaotic 
domain. In order to firmly establish the universality of these results, it is essential to carry 
out numericl studies for a wide variety of interacting particle systems and investigate various 
deformed EGOE's in detail. 
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